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Abstract 

We propose to study the origin of algebraic decay of two-point correlation functions observed 
in glasses, proteins, and quantum dots by their nonlinear response to sequences of ultrafast laser 
pulses. Power-law spectral singularities and temporal relaxation in two-dimensional correlation 
spectroscopy (2DCS) signals are predicted for a continuous time random walk model of stochastic 
spectral jumps in a two level system with a power-law distribution of waiting times ip(t) ~ t~ a ~ l . 
Spectroscopic signatures of stationary ensembles for 1 < a < 2 and aging effects in nonstationary 
ensembles with < a < 1 are identified. 
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Exponential relaxation profiles and correlation functions observed in systems coupled 
to a bath with a short correlation time are signatures of fast memory loss. Such normal 
relaxation can be described by Markovian (e.g. Langevin, Fokker-Planck and master) equa- 
tions of motion; Multipoint Green's functions may then be factorized into products of two 
point functions and carry no additional information. However, correlation functions in many 
systems characterized by a complex free-energy landscape with a broad distribution of dy- 
namical barriers jl], may acquire long stretched-exponential or algebraic tails Q]. Anomalous 
variations of fluorescence emission (blinking) has been observed in biomolecules (3f, poly- 

n n n a 

mers [4J, quantum dots [5|, |6|, 13] and glasses. The fluorescence trajectories are commonly 
analyzed using a continuous time random walk (CTRW) model of spectral diffusion. In this 
letter we demonstrate how spectral lineshapes obtained from two-dimensional correlation 
spectroscopy (2DCS) jg] may be used to probe complex anomalous relaxation processes in 
the condensed phase. 2DCS techniques are femtosecond optical analogues of NMR, that 



have been successful 



y employed towards the study of the structure of peptides [91], chemical 



exchange in liquids 10J and exciton migration in photosynthetic antennae 2DCS pro- 



vides a bulk alternative to single molecule measurements of multipoint correlation functions 
[2I useful for testing microscopic dynamic model. Connection to anomalous relaxation re- 
quires a consistent theory of 2DCS signals with Nonmarkovian spectral fluctuations, which 
is the goal of this letter. 

The CTRW model is defined by a waiting time probability density function (WTDF) 
i/j(t) for stochastic jumps between various states. All memory is erased at the time of 
the jump. This renewal (resetting) property makes it possible to compute all statistical 
measures even in the absence of a Markovian description for the probability distribution, 
and provides a convenient formalism for describing long-term memory effects. When ifj(t) 
has an exponential form ip(t) = exp (—t/K\) (where K\ = tip(t)dt is mean waiting time) 
the system may be described by a Makovian master equation, and the relaxation is normal. 
However, models with a long time algebraic decay ip(t) ~ t _a_1 , show anomalous phenomena 
at long times when the second moment of ip{t) diverges < a < 2. We consider a random 
walk which is observed starting at time 0. The WTDF of the first jump ip'(t) may differ 
from if)(t) since it depends on how the system was prepared before t = 0. It must be 
treated with care, since it strongly affects the nature of the anomalous ensemble. Stationary 
processes must satisfy microscopic reversibility which implies that ip'(t) is given by a product 
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of the survival probability <p(t) = ip(t')dt' that no jump had occurred for time t and the 
equilibrium density of jumps 1/ki resulting in ip'(t) = <p(t)/Ki [l2j. A stationary CTRW 
is thus only possible for 1 < a < 2 where the first moment K\ is finite. The resulting 
power-law relaxation is slow, but eventually the system reaches an equilibrium [ljj. For < 
a < 1, K\ diverges. The system never equilibrates and shows aging effects (i.e. dependence 
on the initial observation time). Aging phenomena are commonly studied by preparing a 
nonstationary ensemble where all particles are assumed to make a jump at the time origin, 
so that ip'{t) = ip{t). Signatures of aging were observed in fluorescence blinking of single 
CdSe quantum dots, analyzed within CTRW and yielded a ~ 0.5 . In contrast, the 

anomalous multipoint correlations observed in fluorescence trace of conformation dynamics 
of flavin proteins showed symmetries due to microscopic reversibility indicative of a 
stationary process. 

We consider a system undergoing anomalous stochastic jumps between two states a and 
b. The system is further coupled to a two-level chromophore with a ground \g) and an 
excited state |e) and dipole moment /x eff , causing its transition frequency Q eg to undergo 
stochastic fluctuations 5Q eg (t). 5Q eg can assume two values Qq and — Qq corresponding to 
the system in state a and b respectively. This is known as the two state jump (TSJ) model of 
spectral lineshapes [15j . Observable quantities are obtained by averaging over the ensemble 
of stochastic paths of 8Q eg (t). 

We propose to probe this complex dynamics through the four-wave-mixing signal gener- 
ated by the coherent response of the chromophore to three short optical pulses with wavevec- 
tors ki, k 2 and k 3 at times 0,t\,ti + 1 2 . The signal detected at time ti + t 2 + 1 3 is described 
by the nonlinear response function which depends on the three experimentally controlled 
parameters ti,t2,t 3 . The photqri^echo signal generated in the ki = — ki + k 2 + k 3 phase- 
matching direction is given by |l( 



^i(t 3 ,t 2 ,ti) = 2(i/h) 3 ^Mh)o(t2)e(t 1 )e in ^-^ 



eg 

x e~ r( * 1+i3) /exp 



rti+t2+ta 

-i / 5{l eq (T)dr 

Jt!+t 2 



exp 



17} 



fti 

5Q eg (r)dT 



(1) 



with 7] = 1. We have added a homogenous dephasing rate T. We shall display this signal by 
2D Ui,u 3 correlation plots 

Sj^h,^) =-Im\ [ ^(hMihY^+^dhdh (2) 
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where t% is held fixed. A different signal, J^ji, in the kn = ki — k 2 + k 3 direction is sim- 
ilarly given by setting rj = — 1. We further considered the combination Sa(&3, t 2 , ^1) = 
Si(uj 3} t 2} —0Ji) + Sjj(lo 3i t 2 , ooi) which shows absorptive peaks [it}]. An optical coherence is 
encoded in the system during t\ and probed during t 3 . By varying the delay t 2 between 
these periods (which is typically much longer than t\ and t 3 ) we can explore correlations be- 
tween dynamical events of the stochastic system. Multipoint correlation functions ordinarily 
obtained in single molecule spectroscopy can thus be obtained from bulk measurements. 

Microscopic reversibility in stationary ensembles implies that S"i(t 3 , t 2 , ti) = ^/(ti, £2, £3) 
and yuits, t 2 , ti) = yn(t 3 ,t 2 ,ti), resulting in the following symmetry of the lineshape 

Sa(u3, h, ^l) = Sa(vi, fo, u 3 ). (3) 

When the system has no memory during the t 3 interval regarding its state during t\ 
Si, and Sjj may be factorized as Si(t 3 ,t 2 ,ti) = 2(i/h)K(t 3 )K*(t 1 ) , and Su{H 1 t2^i) = 
2(i/h)K(t 3 )K(t 1 ). Here K(t) = (i/h)i2 2 eg e^ r+ne9 ^ (exp[-i ft 5fi e9 (r)dr]), is the linear re- 
sponse function. The correlation signal reduces in this case to the product of the linear 
absorption Wa{oj) = Im / °° K(t) exp[icut]dt lineshapes 

hS A (uJ 3 , t 2 -> 00, wx) = AWaMWaM- (4) 

Algebraic memory will result in a slow convergence to this asymptotic lineshapes. In addi- 
tion, the spectra will diverge at certain frequencies where the factorization (Eq. (j3J)) does 
not hold, as will be shown below. 

The response functions will be calculated by introducing a 2 x 2 matrix G v in a, h space 
whose jl element gives the contribution of paths with initial state I and final state j by 
averaging over I and summing over j. 

^u{h, t 2 , h) = {Gu)ji (*s, t 2 , h) \p\ (t = 0) (5) 
jl 

In the case of Markovian relaxation, G u is given by a product of three Green's functions 
representing the time evolution during the t\, t 2 and t 3 intervals. These Green's functions 
can be calculated by solving the stochastic Liouville equation, which combines the Liouville 



equation for coherence evolution with a rate equation for the jumps 18j. The anomalous 
four-point G v , in contrast, may not be factorized in this manner. Calculating it requires 
some bookkeeping of jump events. For each of the three time intervals ti,t 2 , t 3 we distinguish 



between two possibilities; either there was no jump or there was at least one jump during that 
interval. According to this classification, G u is given by a sum of 8 terms G™, (rn — 1, . . . , 8). 
G™ are constructed by convoluting matrix factors for propagation between the first and the 
last jump in the ti interval (if any jump occurred in ti) and factors for the coherent evolution 
between the consecutive jumps in different intervals (last in earlier and first in later) [121 ] . 
The calculation is conveniently made in Laplace domain where convolutions become simple 
multiplications and the integral equation for the propagator factor becomes an algebraic 
equation. Making use of the causality of the response functions (Eq.Q), the 2D lineshapes 
(Eqj2]) were obtained by analytical continuation of the Laplace domain response functions 
calculated in Appendix D of Ref. [3] by setting si = T+irj(uJi— Vt eg ), and s 3 = T—i(uj 3 —Vl e g), 
where Sj is the Laplace variable conjugated to tj. 

We first consider a stationary ensemble with anomalous relaxation where 12] : 



ip(s) is the Laplace transform of ip{t). This WTDF has a mean waiting time K\ and the 
long time t » (kik^ 1 ) 1 ^ asymptotics ipw(t) ~ K^ l Ki/t a+l . In all plots we use the 
dimensionless frequency units (ujj — Q eg )/Qo by setting Q eg = 0,Q = 1. 

No jumps occur during t± and t% in the slow fluctuation QqKi » 1 limit. The linear 
lineshape has two peaks at frequencies to = ±1, and the 2D spectrum for short delay t 2 ~ 
consists of two diagonal peaks (the system is in the same state during t\ and £ 3 ) at (u>i,uj 3 ) 
(1,1) and (-1, -1) centered at these two frequencies. For normal relaxation the peaks are 
Lorentzian. The anomalous model shows the same peak pattern, but the peaks are divergent 
due to long-tailed correlations. For T = the linear absorption peaks diverge as Au) a ~ 2 where 



Auji = iO\ — Vt eg — Vt, 
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In Fig 1 we display the Si, Sn and Sa signals for slow fluctuations at t 2 = 0, and T = 0. 
All panels show prominent two diagonal peaks at (1,1) and (-1,-1). 5*^ is simpler due to 
interference of Si and Sn, which have opposite signs in specific Aui,Au 3 quadrants, and 
substantially cancel. Along the uj\ = ±1 and u>s = ±1 lines, Si and S^/diverge as ~ Au^~ 2 
and ~ Au)i~ 2 respectively. Sa is finite along u)\ — ±1 (except for the u 3 = ±1 peaks), but 
is nondifferentiable with respect to u\. The uj 3 = ±1 line has the same analytical structure. 

The most rapidly divergent term which represents the paths that has no jump during 
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the entire time t± + ti + t 3 is given by (for the (1,1) diagonal peak) 



Q 8{ + r\ \ 4 /4r 
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h 



; Im<j)'(T - iAu 3 ) 



(7) 



+ 2 7— 2^0'(r - iAcUi) 



Q-3 



For T > Eq© is regular, for T = it diverges as S A « ^^T 1 sin [tt(2 - a)/2]A^3 
(ACU3 — > 0,Ao;i = 0). This is illustrated in the linear log-log plot in Fig 2A showing in S A 
vs. co>3 for u)\ — 1. The a — 3 exponent is shown in the inset. 

For fast fluctuations Qqk,\ << 1 (not shown), the state of the system rapidly changes 
during the t 1? t 3 intervals. For normal relaxation, both ID and 2D lineshapes then consist 
of a single Lorentzian peak at the average frequency (motional narrowing). For anomalous 
dynamics the survival probability in the initial state is substantial even for fast fluctuations 
and the two peak divergencies of the slow limit (Fig 1) are still retained. An additional 
central peak at (0, 0) shows up. 

The variation of S A with delay time t% in the slow fluctuation limit is displayed in Fig 
3. We see a buildup of (finite, as shown at right panel) off-diagonal cross-peaks, whose 
lineshapes are dominated by paths where the system is in a different state during t\ and t 3 . 
The (-1,1) peak represents paths where the system was in state a during t± and b during 
t 3 . The (1, —1) peak represents the reverse sequence. Note that this 2DCS equilibrium 
measurement records subensembles of trajectories without having to perturb the system or 
perform a single molecule measurement. Contour elongation of both diagonal and cross 
peaks along the uj\^ = ±1 directions is a signature of the long-time memory and Sa may 
not be factorized as in Eq. (j3J) at these frequencies. Outside these regions, or for finite V the 
asymptotic lineshape Eq. (j4j) is approached algebraically as £2 is increased. By including 
a finite dephasing rate T, peak divergencies are cured and the lineshapes do not differ 
significantly from the normal relaxation case. However, the population evolution during ti 
is not affected by dephasing and the cross-peak dynamics is still anomalous. To study how 
the (1,-1) cross peak grows to its long time value (Eq. ([3])) we have considered the quantity 
Z{t 2 ) = \S A (wx,t 2 ,u 3 ) —S A (wi, 00, u 3 )\/\S A (wi, 0,073) -S A {ui, 00,073)! for u x = 1, u 3 = -1. 
Straight lines in log-log plots of Z(t%) vs. £2 (Fig 2B) indicate algebraic relaxation. The 
asymptotic exponent (slope in Fig 2B) approaches a — 1, which is identical to that of the 
frequency correlation function (5Q eg (t)5Q eg (0)) [l^ . 
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Finally we discuss aging effects for < a < 1 case, where k% diverges. This describes 
nonergodic nonstationary processes which never equilibrate, and their properties change with 
time [19J. We model it by a random walk which is started by a jump at some fixed time. The 
response function depends on the time elapsed from the start of a random walk t in to the 
first laser pulse at [20j, l2l| . This effect is fully described by allowing ip'(t;t in ) to depend 



on t in . In Laplace space we find ip'(s;s in ) = [ip{s in ) - ^(s)]^(s in )/[l - ^(s in )](s - s 



in I 



22( | . The nonlinear lineshapes measured in the time-domain provide a direct measure of 



response function, unlike frequency- domain measurements which involve stationarity in order 
to connect with response function [15j]. However, due to the lack of equilibration, averaging 
over consecutive pulse sequences may depend on the repetition rate: A fresh ensemble 
with the same ip'(t) must be prepared before each application of the pulse sequence for the 
response function to be physically meaningful. 

We have calculated the response functions for the model ^jv(s) = 1/[1 + (ks)"], where 
ip(t) ~ (/t/t) 1+a and < a < 1. In Fig IB the Sa lineshape is shown for £j n = 0, i.e. 
the random walk is started each time the first pulse interacts with the sample and ip'(t) = 
ip(t) [19| . Microscopic reversibility breaks down for nonstationary processes and obviously 
the symmetry (Eq. Q) is violated. The nonstationary anomalous ensembles starts with 
some jump rate which depends on the initial preparation. At long times it approaches the 
equilibrium value which vanishes for the present model 2l|. The higher rate during 

the (earlier) t\ interval compared to results in broader peaks along uj\ axis (compared to 
003). Contour elongation along the uji^ = ±1 axes is again observed with divergent diagonal 
peaks of complex structure. Our simulations demonstrate that two-dimensional correlation 
plots of the signals obtained from the response of the system to sequences of multiple laser 
pulses carry specific and direct signatures of complex dynamics. 

The support of the MSMT CR (MSM 0021620835), NSF (CHE-0446555) and NIH 
(2RO1GM59230-05) is gratefully acknowledged. 
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Figure captions 



Fig 1 (Color Online) (A): The 2D Sj(—Ui,u! 3 ) (left), Sh(l)i, u> 3 ) (middle) and Sa(ui,u> 3 ) 
(right) signal (Eq.©) for the WTDF (Eq.flgJ) at t 2 = 0,for a = 1.2, and ««/«! = 0.25, 
fio^i = 2.0, fio = 1; &>eg = 0. (B): The 2D Sa signal for the nonstationary random 
walk a = 0.5, t 2 = 0, t in = kQ = 2.0, Q = 1, f2 e9 = 0. 

Fig 2 (Color online) Panel (A): Peak divergence along Aui = 0. t 2 = 0, a = 1.2(dotted), 
1.5 (dashed), and 1.8 (solid). Other parameters are the same as in Fig 1. Inset: 
peak exponent 7 = ^°g[^(A^i|=o^=o,A^)] _ Pand ( B ). z ^ for the WTDF ( Eq . 

K a / K% = 0.5 , Q Ki = 1.0, a = 1.2(dotted), a = 1.5 (dashed), and a = 1.8 (solid), 
flo = 1. Power-law growth of cross peaks is seen. 

Fig 3 (Color Online) Variation of the 2D S A signal (Eq.©) for the WTDF (Eq. ©) with 
the time delays (left to right) t 2 = 0, 2k±, 10%. Other parameters same as Fig (1A). 
Right panel: Cross peak for t 2 = lO^i on an expanded scale. 
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Fig 1 
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Fig 3 
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